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ABSTRACT 

The covering factor of Compton-thick obscuring material associated with the torus in active galactic 
nuclei (AGN) is at present best understood through the fraction of sources exhibiting Compton-thick 
absorption along the line of sight (Ah> 1.5 x 10 24 cm -2 ) in the X-ray band, which reveals the average 
covering factor. Determining this Compton-thick fraction is difficult however, due to the extreme 
obscuration. With its spectral coverage at hard X-rays (>10 keV), NuSTAR is sensitive to the AGN 
covering factor since Compton scattering of X-rays off optically thick material dominates at these 
energies. We present a spectral analysis of 10 AGN observed with NuSTAR where the obscuring 
medium is optically thick to Compton scattering, so called Compton-thick (CT) AGN. We use the 
torus models of Brightman & Nandra which predict the X-ray spectrum from reprocessing in a torus 
and include the torus opening angle as a free parameter and aim to determine the covering factor 
of the Compton-thick gas in these sources individually. Across the sample we find mild to heavy 
Compton-thick columns, with N-r measured from 10 24 10 26 cm -2 , and a wide range of covering 
factors, where individual measurements range from 0.2-0.9. We find that the covering factor, / c , 
is a strongly decreasing function of the intrinsic 2-10 keV luminosity, Lx, where / c = (—0.41 ± 
0.13)logio(Lx/ergs -1 )+18.31 ± 5.33, across more than two orders of magnitude in Lx (10 415 10 44 
ergs -1 ). The covering factors measured here agree well with the obscured fraction as a function of 
Lx as determined by studies of local AGN with Lx> 10 42 ' 5 ergs -1 . 

Subject headings: galaxies: active — galaxies: individual (NGC 424, NGC 1068, 2MFGC 2280, NGC 
1320, NGC 1386, NGC 3079, IC 2560, Mrk 34, NGC 4945, Circinus) — galaxies: 
nuclei — galaxies: Seyfert — X-rays: galaxies 


1. INTRODUCTION 

The average covering factor of the obscurer in active 
galactic nuclei (AGN), understood to be a torus-like 
structure, is represented by the ratio of obscured to un¬ 
obscured AGN, or the obscured fraction. The obscured 


A 

Lawrence 

1991 Simpson 

20051, the 

Lusso et al. 

2013; Gu 2013 

Asset et ai. 

-rays (e.g. 

Lisaliti et al. 

IS 

99) and is 


largely understood to be dependent on the power of the 
central source (e.g. Ueda et al. 2003 La Franca et al. 


2005 

|Akylas et al. 

2006; Tueiler et al. 2()08| | Basinger 

2008 

Beckmann et al. ||2009; 

Akylas & Georgantopoulos 

2009 

Burlon et ai. 12011 Brig 

ltman & JN andra 2011b Va- 


sudevan 
that this may 


be an observational effe ct (e.g. Lawrence 
fe Elvis|[2010 Mayo fe Lawrence||2013 1. 


The fraction of AGJN where the obscuration is so ex- 


























































2 


Brightman et al. 


treme that it is optically thick to Compton scatter¬ 
ing {Nh> 1.5 x 10 24 cm -2 ), so call ed Compton -thick 


(CT) AGN, is less well know n (e.g. Burlon et al. 2011 
Brightman & Nandra 2011b), with onl y a tew tens o 


bonahde CT AGIN know n locally (e.g. Goulding et al. 
2012| |Gandhi et al.||2014[). The dependency ot the cov- 
ermg factor ot Compton-thick gas on the power of the 
AGN is relatively unknown, even in the local universe, 
hindered by low number statistics. At the higher red- 
shifts probed by deep extragalactic surveys by Chandra 
and XMM-Newton, the redshifting of the Compton-hump 
into the bandpasses o f these telescopes aid s identification 
of these sources (e.g. Comastri et al.pOlT). In addition, 


the large volumes probed yield ot order hundreds to thou¬ 
sands of sources overall such th at significant numbers of 


CT AGN can be uncovered (e.g. Brightman fe Ueda 2012 
Buchner et al. 2014 Brightman et al. 2014). INeverthe- 


less, the low count nature ot these sources means that 
constraints on the Nh and Lx are still relatively poor. 

The ability to determine covering factors for individual 
sources is needed to solve many of the issues mentioned 
above and is needed in order to carry out a detailed study 
of what physically affects the covering factor. Estimating 
the covering factor in individual sources can be done by 
determining the ratio of emission reprocessed in the ob¬ 
scuring medium to the intrinsic emission, which is best 
done in unobscured sources (i.e. those viewed with a 
clear line of sight to the nucl eus) where the intrinsic emis- 


sion is directly visible (e.g. iTreister et al. 

2008 

Gandhi 

et al. 

2009 

Lusso et al. 

2013[Gu 2013). 

However the 


reprocessed emission, making them difficult to disentan¬ 
gle and furthermore disk reflection is difficult to separate 
from distant reflection. In obscured sources (i.e. those 
viewed through thick material) the reprocessed emission 
dominates, though estimates of the intrinsic emission are 
challenging due to the obscuration itself. 

X-ray spectral analysis extending to high energies is 
ideal for determining the covering factor due to the fact 
that X-rays above ^3 keV penetrate all but the most 
extreme obscuration, allowing a good estimate of the in¬ 
trinsic power. Furthermore, in CT AGN where the ob¬ 
scuring medium is optically thick to Compton scattering, 
the scattering of X-rays within the medium can reveal 
the covering factor of the gas for such sources, especially 
evident above 10 keV. This has been done for a hand- 
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the AGN covering factor has yet to be carried out. 

The recently launch ed Nuclear Spectro scopic Telescope 
Array ( NuSTAR , Harrison et ah 2013) is sensitive in 
the 3-79 keV band. Its significantly improved sensitivity 
above 10 keV with respect to previous telescopes makes 
it the ideal instrument to measure the strength and shape 
of the Compton reflection hump and thus determine the 
covering factor in a large number of AGN. As part of its 
baseline mission, NuSTAR has observed ~100 AGN from 
the hard X-ray selected Swift /BAT survey for ~20 ks 
each, as well as longer observations of several well known 
obscured AGN such as NGC 1068 and NGC 4945. This 
sample, which has also been studied at many other wave¬ 
lengths, is ideal for investigating the covering factor and 


how it varies. In this paper we present an initial analysis 
of a small sample of 10 of these sources to investigate 
how well the covering factor can be determined from X- 
ray spectra. While this sample is small, it was selected 
in order to cover a wide range in Lx, paying particular 
attention to low (L x ~ 10 42 ergs -1 ) and high (L x ~ 10 44 
ergs -1 ) luminosity sources to make it as representative 
as possible. In Section 2, we introduce X-ray spectral 
torus models and make comparisons between them. In 
Section 3, we describe the sample and the data analysis. 
In Section 4, we present our spectral fitting results. In 
Section 5, we discuss potential biases and systematics. 
In Section 6, we compare our results to previous results. 
In Section 7, we explore what physically influences the 
AGN covering factor. Finally, in Section 8, we present 
our conclusions. We assume a flat cosmological model 
with ffo=70 km s -1 Mpc -1 and Ga=0.73. For mea¬ 
surement uncertainties on our spectral fit parameters we 
present the 90% confidence limits given two interesting 
parameters (Ay=4.61). 

2. X-RAY SPECTRAL TORUS MODELS 

We first briefly summarize the development of X-ray 
spectral torus models, detailing their similarities and dif¬ 
ferences. We then illustrate the differences in the X-ray 
spectra resulting from these models by fitting simulated 
spectra created from one model to other sets of public 
models and comparing the results. 


2.1. Model details 

In recent years a suite of new X-ray spectral models 
that describe the reprocessing o f X-rays in a torus-shaped 
medium have been published fllkeda et al. 2009 Mur- 
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A 20141). These models employ Monte-Carlo techniques 


to calculate spectra, simulating photoelectric absorption, 
fluorescence an d Co mpton s catte rin g and buil d on pre¬ 
vious work by Matt et al. J 1991 b_ Leahy fe Creighton 
dl99 3|), iGhisellinf et al.| (|7ffij4[) JAandra &; Georgel ( 1994[) 
and Yaqoob 1 1997|). The aforementioned models differ 
mostly in the geometry of the torus considered and the 
treatment of the different components, be it direct, scat¬ 
tered or line emission. The underlying physics of the 
models is th e same, however, w ith photoelectric c ross¬ 


sections from Verner et al. (1996), abundances from Am 
ders fe Gre vcsse (1989) and fluorescent yields fromlBam- 
bynek et af. 11972) ~However, the physical geometry and 


treatment of Com pton scattering d iffers somewhat. The 
latest edition from Liu & Li (20141 is the first to explic¬ 


itly simulate a cl umpy torus. 


The model of Ikeda et al. (2009) considers a spher¬ 
ical torus, which consists of a spherical distribution of 
matter with a biconical void. The column density, , 
through the torus is a function of the inclination angle, 
where the Nr is maximum when the torus is seen edge- 
on (90° inclination angle). This decreases to zero as the 
inclination angle decreases and reaches the opening an¬ 
gle, at which point the source can be seen unobscured. 
The opening angle is a free parameter with a range of 
0-70°. The model separates into three components: the 
direct zeroth-order transmitted component, the absorbed 
reflected component and the unabsorbed reflected com¬ 
ponent. 
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Murphy & Yaqoob (2009) consider a torus with a circu¬ 
lar cross section and a fixed opening angle of 60° (cover¬ 
ing factor of 0.5). This model, known as MYTORUS, also 
has separate components for the direct, scattered and 
line emission, motivated so that time variability between 
these different emissions can be studied. The line com¬ 
ponents are separated so as to account for instrumental 
system atic effects on the line energies. Both MYTORUS 
and Ikeda et al. (2009) consider Ah up to 10 25 cm -2 . 

A new 'imple mentation of MYTORUS was introduced by 
Yaqoob (2012) in which the various model components 
can be used independently in so called ‘decoupled’ mode, 
giving it added flexibility. In this way, it can be used to 
model a clumpy distribution of matter with an arbitrary 
covering factor. 


The model of Brightman & Nandra (2011a) als o as¬ 
sumes a sph erical torus with a biconical void, as in |Ikeda| 
et al. (2009). However, the line of sight Ah through the 


torus is constant and not dependent on the inclination 
angle. This model also allows for a variable opening an¬ 
gle, ranging from 26-84°, and has the added advantage 
of extending up to Ah=10 26 cm -2 , thus allowing inves¬ 
tigations of the most extreme obscuration. We refer to 
this model simply as t h e TORUS model. In addition to 
this TORUS model, Brightman & Nandra (2011a) present 
a model for the special case where the source is fully 
covered (a spherical geometry) and includes variable el¬ 
emental and iron abundances with respect to hydrogen. 
We refer to this model as the sphere model. 


Recently Liu & Li (2014) have investigated the effect 
of a clumpy medium, distributed in a toroidal geometry, 
on the emergent X-ray spectra. The parameters they 
consider are the volume filling factor of the clumps, the 
number of clumps along the line of sight, and the line of 
sight Ah- They also only consider Ah up to 10 25 cm -2 . 

In our in vestigation of the cover ing facto r in CT AGN, 
we use the Brightman & Nandra (2011a) TORUS model 
because it includes a variable torus opening angle and 
extends to a higher Ah range. Henceforth we will refer 
to the half-opening angle of the torus, measured from 
the polar axis to the edge of the torus, as 9 tor . The 
covering factor of the torus is simply calculated from 9 toI 
following / C =cos(0 t or)- Figure [l] shows spectra from this 
model, where the left panel shows a range of 9 tor values 
(equivalent to a range of f c values), and the right panel 
shows a range in inclination angles (0; nc )• The effect of 
the variation of the opening angle can clearly be seen in 
the left panel over the entire 3-79 keV spectral range of 
NuSTAR , in particular in the strength of the Compton 
hump at ~ 20 keV. This illustrates how NuSTAR is well 
suited for studying the covering factor in CT AGN. 

2.2. Model comparisons 

As these X-ray torus models are still fairly novel, inter¬ 
preting their results when applied to real data is not yet 
fully tested or understood. In order to interpret our re¬ 
sults using the TORUS model and to make inferences from 
them, we aim to better discern how this model relates 
to other commonly used models. We therefo re compare 
the TORUS model to the PEXRAV model of Magdziarz 
& Zdziarski (1995), which describes the Compton re¬ 
flection ofX-rays off a semi-infinite slab of cold mate¬ 
rial. The PEXRAV model has been extensively used in 
the literature for describing the effect of Compton scat¬ 


tering/reflection in AGN spectra and for fitting the spec¬ 
tra of reflection-dominated, CT AGN. Furthermore, we 
compare the TORUS model to MYTORUS, described above. 
These two models have differing geometries and thus a 
comparison highlights the effects of geometry on the X- 
ray spectra. MYTORUS has been used frequently in the 
anal ysis of NuSTAR spec t ra of heavily obscured AGN 
(e.g. iPuccetti et al. 2014||Gkuwihijtjil. 2014 JA revalo 


ILE 

uab 


et al.|2()TT ' Balokovic et al. 201 1: iBauer et aC poTH i , and 
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parison to the 
torus models o: 
public. 

To make the comparisons, we simulate NuSTAR spec¬ 
tra from the TORUS model for Ah=10 24 , 10 25 and 10 26 
cm -2 and for 6^=40°, 60° and 80° for F = 1.9 and 
a 3-79 keV flux of ~ 10 -11 ergcm _2 s _1 . We simulate 
NuSTAR. spectra with 100 ks exposures without count¬ 
ing statistics in order to only consider the effect of the 
model parameters on the spectrum, and we use response 
and background files from a randomly selected observa¬ 
tion, that of NGC 1320 (see Tablejl] obsid 60061036004). 
As a consistency check, we fit the TORUS model to the 
simulated data, where we find that the input parameters 
are recovered. 

We subsequently fit each simulated spectrum between 
3-79 keV with the pure reflection component of PEXRAV 
with the high-energy cut off in the continuum set to the 
maximum (l x 10 6 eV), the abundances at solar and 
the cosine of the inclination angle at 0.45. We add a 
narrow Gaussian component at 6.4 keV to model the 
Fe Ka line, since unlike TORUS, PEXRAV does not self- 
consistently include Fe fluorescence. For Ah=10 24 cm -2 , 
an additional absorbed power-law component is required, 
which we model with ZWABS x POWERLAW, the pho ton 
index of which is fixed to that of PEXRAV. Figure |2.2| 
shows how the best fit model compares to the simulated 
spectra. For TORUS Ah= 10 24 cm -2 , where the optical 
depth to Compton scattering is just below unity, the fit¬ 
ted models agree with the TORUS model, with no obvious 
deviations from unity in the data-to-model ratio. This is 
also the case for all Ah values with a 60° opening angle. 
However, at small (40°) and large (80°) opening angles, 
when fitting with PEXRAV there are obvious deviations in 
the data-to-model ratio above 10 keV. This implies that 
PEXRAV is not able to reproduce the varying high-energy 
spectral shapes produced by a torus geometry. 

Following this we fit MYTORUS to these simulated spec¬ 
tra in both coupled and decoupled mode. In coupled 
mode the parameters of the scattered, direct and line 
components are fixed to one another and thus the cover¬ 
ing factor is 0.5. For the decoupled mode we combine the 
direct absorbed component with two scattered and line 
components, one where the inclination is fixed to 0°, the 
other where it is fixed to 90°, with all parameters linked 
with the exemption of the normalisations. This imple¬ 
mentation represents backwards scattering and f orwar ds 
scatt ering respectively, which as described by Yaqoob| 
(2012), mimics an absorber with a clumpy distribution 
of matter with a variable covering factor. The data-to- 
model ratios for the fit with MYTORUS in decoupled mode 
are shown in Fig. [3] These show that the flexibility of 
MYTORUS in this mode can reproduce the spectral shape 
of the TORUS model for a wide range of parameters, show- 
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Figure 1 . X-ray TORUS spectra for N^=2 X 10 24 cm -2 , T=2 and different torus opening angles, 9toi (left), where the equivalent covering 
factor is in parenthesis, and different inclination angles, f! mc (right). The variation in spectral shape with 0tor above 10 keV is clear, 
illustrating the capability of NuSTAR to determine the torus opening angle for CT AGN. For a given torus opening angle, changes in the 
inclination angle primarily affect the < 10 keV X-ray spectrum. 


ing deviations only for the extreme case of ^tor^O 0 and 
1 Vh= 10 25 cm -2 and above. 

Since the power-law index, T, is a good indicator of the 
overall shape of the spectrum, we compare that parame¬ 
ter for the fitted models to the value of T adopted by the 
input model. We do not aim to test if the fitted models 
can recover the intrinsic T, since the models have differ¬ 
ent geometries, rather we explore the differences in T in 
order to gain insight into the various spectral shapes pro¬ 
duced by each model. This will also allow better inter¬ 
pretation of the parameters obtained when fitting these 
models to real data. Figure [4] shows how the best-fit T 
from PEXRAV and MYTORUS compare to the input value 
of 1.9 from the TORUS model given the input Ah and 
9 tm parameters. For all models, the measured T is con¬ 
sistent with the input T at Ah=10 24 cm -2 . For PEXRAV, 
this diverges to lower T values at higher Ayr for 40° and 
60° opening angles, indicating that the TORUS model is 
producing a stronger Compton hump than PEXRAV for 
these parameters. MYTORUS in coupled mode is roughly 
consistent at 40°, whereas 60° shows no difference. This 
is expected since the opening angle in this model is in¬ 
deed 60°. For 80°, both PEXRAV and MYTORUS mea¬ 
sure rather larger T values, indicating that the Compton 
hump is weaker in the TORUS model for large 9 tOT - 

For MYTORUS in decoupled mode, as in coupled mode, 
the measured T from the simulated spectra agree with 
the input value within the measurement uncertainties for 
the case of 40° and 60°. The measurement uncertainties 
are higher in decoupled mode due to the larger number 
of degrees of freedom with this implementation. For the 
case of 80° there is a large discrepancy between the input 
and recovered F values. 

As described above, the decoupled implementation of 


MY TORUS al l ows fo r an arbitrary covering factor. Follow¬ 


ing Yaqoob (20121, Puccetti et al. (2014) estimated the 


covering factor for INGC 4945 from / c ~ 0.5 x Asgo/Azgo, 
the ratio of the normalizations of the direct and scattered 
components at 90° inclination angles. We can directly 
compare this estimation with our input covering factors. 
Our input 9 t0 i values correspond to f c = 0.77, 0.50 and 
0.17. From the above formulation, we find the cover¬ 
ing factors estimated by the decoupled implementation 
of MYTORUS to be 0.41, 0.09 and 0 respectively, which 
are much lower than the input values. It does however 
recover the input trend of decreasing covering factors. 
The decoupled implementation of MYTORUS only mim¬ 
ics a free covering factor and thus an agreement was not 
necessarily expected. 

Along with th e PE XRAV deviations in spectral shape 
shown in Figure [2~2| the above analysis has shown that 
PEXRAV will also systematically obtain different spectral 
parameters compared to the TORUS model. We have also 
shown that the opening angle is an important parame¬ 
ter in determining the spectral shape, and thus models 
with fixed opening angles, such as MYTORUS, will also 
systematically obtain different parameters. We therefore 
conclude that slab reflection models should not be used 
for fitting th e high-energy X-ray spect ra of CTAGN, as 
concluded by Murphy & Yaqoob (2009), and that ideally 
spectral models with the covering factor of the Compton- 
thick gas as a free parameter should be used. 

3. SAMPLE PROPERTIES AND DATA ANALYSIS 


Our sample consists of 10 local (z < 0.1) CT AGN ob¬ 
served by NuSTAR: NGC 424, NGC 1068, 2MFGC 2280, 
NGC 1320, NGC 1386, NGC 3079, IC 2560, Mrk 34, 
NGC 4945 and the Circinus galaxy. With the exception 
of 2MFGC 2280, these sources were selected to be known 
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Figure 2. Data-to-model ratios of simulated NuSTAR spectra using the TORUS model (red data points) fitted by PEXRAV plus an absorbed 
power-law for A r H=10 24 cm -2 , for three TORUS JVh values (JVh = 10 24 , 10 25 and 10 26 cm -2 ) and three values of $tor (40°, 60° and 80°). 
For TORUS Nn= 10 24 cm -2 , where the optical depth to Compton scattering is just below unity, the fitted models agree with the TORUS 
model. This is also the case for all IVh values with a 60° opening angle. However at small (40°) and large (80°) opening angles, PEXRAV 
cannot reproduce the shape of the TORUS model. 
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Figure 3. Same as Fig. |2.2[ but simulated spectra fitted by MYTORUS in decoupled mode. 
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Figure 4. Comparison of the power-law index, F. determined 
from fitting PEXRAV and MYTORUS to simulated NuSTAR spectra 
using the TORUS model (r=1.9, dashed line) for 7Vh=10 24 , 10 25 
and 10 26 cm -2 and three opening angles (40°, 60° and 80°). 


CT AGN and to cover a wide range in X-ray luminosity. 
Particular attention was paid to selecting sources in the 
low (Lx~ 10 42 ergs -1 ) and high (~ 10 44 ergs -1 ) lumi¬ 
nosity regimes, thus enveloping the bulk of the local CT 
AGN population, which is at moderate luminosities, so 
that our sample could be as representative as possible, 
despite its small size. CT AGN are selected since the 
obscuration is required to be optically thick to Compton 
scattering in order for the models to be effective at distin¬ 
guishing the covering factor. The models are degenerate 
with covering factor for Compton-thin media. Previously 
published NuSTAR observations were included in order 
to facilitate comparison among models. 

Observati ons of NGC 424, NGC 1320 and IC 2560 are 
presented in Balokovic et al. (20141 with modelling using 
PEXRAV and MYTORU S . Res ults on NGC 1068 are pre¬ 
sented in Bauer et al. (2014), also with modelli ng using 
PEXRAV and MYTORUS. IVlrlc 34 is presented in [Gandhi] 


et al. (2014) with modelling using MYTORUS and TO RUS. 
IN GC 4945 is presented in Puccetti e t al. (2014) and 
Circinus is presented in Arevalo et ai. (|2014); both are 
modelled using MYTORUS. The NuS'l'AR observations of 
NGC 1386, NGC 3079 and 2MFGC 2280 are presented 
here for the first time. No detailed X-ray spectral mod¬ 
elling has previously been published on 2MFGC 2280, 
and it was selected for study here based on its flat spec¬ 
tral shape in the NuSTAR band, which resembles that of 
the other obscured AGN considered in this work. 


Table □ summarizes the basic observational data for 
this sample. The NuSTAR observations were perfo r med 


for se veral different scientific purposes (see Harrison et al. 
2013 for a description of the NuSTAR survey^ programs) 


and therefore differ significantly in exposure time and 
signal-to-noise ratio. We use the first three existing NuS¬ 
TAR observations of NGC 1068, both observations of 
NGC 1320 and both observations of IC 2560 as these 
sources are not known to be variable in the hard X-rays. 
NGC 4945 is known to be variable. We mitigate the 
effect of the variability by analysing only one NuSTAR 
obsID from this source, of which we choose the first of 
three. We investigate the effect on our conclusions of this 
choice later in the paper. For a ti me resolved-analysis o f 
this source, we refer the reader to Puccetti et al. (2014). 

For NGC 424, NGC 1320, NGC 1386, NGC“3U79, 
IC 2560 and Mrk 34, we use additional archival data from 
XMM-Newton from 0.5-10 keV to provide additional 
spectral constraints. XMM-Newton data were chosen 
specifically over other instruments as XMM-Newton and 
NuSTAR are well matched in terms of effective area and 
spatial resolution in the energy range where they com¬ 
monly observe (i.e. 3-10 keV). No good quality data exist 
for 2MFGC 2280 below 10 keV. NGC 1068, NGC 4945 
and Circinus are part of a different, non-snapshot NuS¬ 
TAR program and hence have longer exposure times and 
higher signal-to-noise spectra and are thus of sufficient 
quality to constrain the torus parameters well without 
additional soft X-ray data. 

The raw data were reduced using the NuSTARDAS 
software package version 1.2.1. The events were cleaned 
and filtered using the nupipeline script with standard 
parameter^] The nuproducts task was used to generate 
the spectra and the corresponding response files. Spec¬ 
tra were extracted from circular apertures centred on the 
peak of the point source emission, with radii between 
30"and 90"(larger radii for higher signal-to-background 
ratio). The background spectra were extracted from re¬ 
gions encompassing the same detector as the source, ex¬ 
cluding the source extraction region and avoiding the 
wings of the PSF as much as possible. Data from both fo¬ 
cal plane modules (FPMA and FPMB) and from multiple 
observations are used for simultaneous fitting, without 
co-adding. Additional details of the data reduction for 
specific sources can be found in the papers listed above. 
The observations of NGC 1386 and 2MFGC 2280, which 
are not p ublished elsewhere, ar e reduced in a fashion 
similar to Balokovic et al. (2014). 

The XMM-Newton data reduction for the observations 
of NGC 424, NGC 1320, NGC 1386 and NGC 3079 were 


described in Brightman & Nandra (2011a), where only 
EPIC-pn data were used, with source events extracted 
from 35" radius circular regions. The events were then 
filtered for background flares when the level of the back¬ 
ground count rate was determined to be twice the level at 
which the excess variance was determined to be zero. The 
XMM-Newton data reduction for the observation of Mrk 


34 is described in Gandhi et al. (20141) and the XMM- 
Newton dat a reduction for 1C 2560 is described in ITilakl 
eFaL] (2008). 


1 The NuSTAR Data Analysis Software 
http: //heasarc.gsfc.nasa.gov/docs/nustar/analysis/ 
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Table 1 

Summary of the observational data used in this analysis. 


Source name 

(i) 

RA 

(deg) 

(2) 

Dec 

(deg) 

(3) 

telescope 

(4) 

obsid 

(5) 

date 

(6) 

exposure 

(ks) 

(7) 

net count rate 
(counts s _1 ) 
(8) 

NGC 424 

17.86511 

-38.08347 

NuSTAR 

60061007002 

2013-01-26 

15.5 

0.052/0.053 




XMM-Newton 

2942301 

2001-12-10 

4.5 

0.236 

NGC 1068 

40.66963 

-0.01328 

NuSTAR 

60002030002 

2012-12-18 

57.8 

0.112/0.107 





60002030004 

2012-12-20 

48.6 

0.111/0.111 





60002030006 

2012-12-21 

19.5 

0.115/0.113 

2MFGC 2280 

42.67746 

54.70489 

NuSTAR 

60061030002 

2013-02-16 

15.9 

0.019/0.018 

NGC 1320 

51.20292 

-3.04228 

NuSTAR 

60061036002 

2012-10-25 

14.5 

0.028/0.023 





60061036004 

2013-02-10 

28.0 

0.025/0.022 




XMM-Newton 

405240201 

2006-08-06 

12.5 

0.081 

NGC 1386 

54.19242 

-35.99941 

NuSTAR 

60001063002 

2013-07-19 

21.2 

0.010/0.011 




XMM-Newton 

140950201 

2002-12-29 

13.9 

0.153 

NGC 3079 

150.49085 

55.67979 

NuSTAR 

60061097002 

2013-11-12 

21.5 

0.062/0.063 




XMM-Newton 

110930201 

2001-04-13 

13.6 

0.092 

IC 2560 

154.07799 

-33.56379 

NuSTAR 

50001039002 

2013-01-28 

23.4 

0.012/0.012 





50001039004 

2014-07-16 

49.6 

0.012/0.012 




XMM-Newton 

0203890101 

2003-12-26 

80.5 

0.006/0.007 

Mrk 34 

158.53580 

60.03111 

NuSTAR 

60001134002 

2013-09-19 

23.9 

0.007/0.007 




XMM-Newton 

0306050701 

2005-04-04 

8.8/22.9 

0.070/0.013 

NGC 4945 

196.36449 

-49.46821 

NuSTAR 

60002051002 

2013-02-10 

45.2 

0.263/0.249 

Circinus 

213.29146 

-65.33922 

NuSTAR 

60002039002 

2013-01-25 

53.8 

1.046/0.984 


Note. — Column (1) gives the source name, Columns (2) and (3) list the J2000 position given in NED in degrees, column (4) gives the 
telescope name, column (5) lists the observation ID, column (6) gives the start date of the observation, column (7) gives the exposure time 
in ks and column (8) gives the net count rate for each instrument, be it FPMA/B for NuSTAR (3-79 keV) or pn/MOS for XMM-Newton 
(0.5-10 keV). Only pn data are used for NGC 424, NGC 1320, NGC 1386 and NGC 3079 and only MOS data are used for IC 2560. 


All spectra were grouped with a minimum of 20 counts 
per bin using the HEASARC tool grppha. We use XSPEC 
version 12.6.0 to carry out X-ray spectral fitting, and 
the x 2 statistic as the fit statistic, with the background 
subtracted. For NuSTAR data, only energies from 3-79 
keV were considered, as the calibration at lower ener¬ 
gies is uncertain, and the NuSTAR response cuts off at 
79 keV due to an absorption edge in the optics. XMM- 
Newton data were used from 1-10 keV, allowing consid¬ 
erable overlap with NuSTAR. 

The primary goal of our analysis is to determine the 
covering factor of the obscuring material in these CT 
AGN using the TORUS model. The AGN emission which 
is described by this model is by far the dominant emission 
in the NuSTAR band for these sources, thus all spectra 
were fitted with the TORUS model as a baseline. We fix 
the inclination angle of the torus to the edge-on position 
of 87°, the maximum inclination angle allowed by the 
model. As described above, the X-ray spectrum above 10 
keV is largely insensitive to the inclination angle, when 
the inclination angle is greater than 0 tor . Fixing the incli¬ 
nation angle thus reduces the number of free parameters 
in the fit and the edge-on choice allows the exploration 
of the full range of opening angles as opening angles can 
only be determined up to the inclination a ngle, after 
which , the source becomes unobscured. | Arevalo et al.| 
(2014) find that when fitting the TORUS model to the 
spectrum of Circinus, that both 9 t0 r and the inclination 
angle can be constrained and provide a better fit to the 
data than if the inclination angle were fixed. We thus 
allow the inclination angle to be free when fitting the 
spectrum of Circinus. 

While the torus emission is the dominant component 
in the NuSTAR band, emission not directly associated 


with the intrinsic AGN emission or the reprocessing in 
the torus, such as soft emission fro m photo-i o nised rnate- 
rial, radiative recombination (Guainazzi & Bianchij2007) 
and Thompson-scattered AGN fight, are common in ob- 
scured AGN and are non-negligible even in the NuSTAR 
band. Furthermore, due to the size of NuSTAR's PSF, 
X-ray sources in the host galax y of the AGN may also 
contribute. For example, Puccetti et al. (2014) finds that 
~ 60% of emission in the 4 — 6 keV band comes from 
extra-nuclear sources within the NuSTAR extraction re¬ 


gion, while Bauer et al. (20141 finds that 28% of the Fe 
Ka emission from N GC 1068 comes from the host galaxy. 
It is expected, however, that these extra-nuclear sources 
contribute far less at higher energies. We therefore add a 
power-law model and a thermal plasma component mod¬ 
elled by APEC to broadly account for these soft excess 
components and extra-nuclear sources. We allow the 
secondary power-law index to vary freely, likewise for 
the temperature of the APEC model. For 2MFGC 2280 
where only NuSTAR data are available, the soft com¬ 
ponents were not statistically required, so they were re¬ 
moved from the fit. No secondary power-law was re¬ 
quired in the fit for IC 2560. For NGC 1068, NGC 4945 
and Circinus where only NuSTAR data were used, the 
APEC model was not required. This is expected since 
this model generally describes emission outside the NuS¬ 
TAR bandpass. The spectra of NGC 424 and NGC 1320 
also did not require the APEC model. 

NGC 1068 has a well known Fe complex at 6-7 keV, 
consisting of neutral Fe Ka and /3 emission, plus ionised 
emission at 6.7 and 6.96 keV. While the TORUS model 
accounts for the neutral emission, we add Gaussian com¬ 
ponents fixed at the energies of the ionised emission to 
account for these lines, where the widths are fixed at 
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sma ll values (1 eV ) . As f ound by Sambruna et al. 


Molendi et al. (2003) for Circinus and Yaqoob 


( 2001 ) 

( 2012 ) 


and 

for NGC 4945, we also iind that the NuSTAR spectra of 
these AGN require an additional Gaussian component 
each to model lines in the Fe complex, fixed at 6.7 keV 
in NGC 4945 and allowed to vary around 6.6 keV for 
Circinus, as required by the data. Additional Gaussian 
components at the energies of the Fe complex were re¬ 
quired to fit the spectra of NGC 1320, NGC 1386 and 
IC 2560. We summarize all the models used for each 
source in Table [2] 

We allow the cross-normalization between the two 
NuSTAR focal plane modules (FPMs) to vary in the fit¬ 
ting to account for instrumental cross-calibration. We 
also allow the cross-normalization between the FPMs 
and XMM-Newton to vary to account for instrumental 
cross-calibration and the fact that the the XMM-Newton 
observations are not simultaneous with NuSTAR. We fix 
the other parameters to each other since there is a known 
agreement between spectral parame ters determ ined from 
XMM-Newton and NuSTAR (e.g. Walton et al. 12013 
2014|. We present the cross-normalizations in Table [2 
The" observations are mainly consistent with the cross¬ 
normalization between the FPMs being unity, with the 
exceptions being NGC 1068, NGC 1320 and Circinus. 
NGC 1068 and Circinus both show that the normaliza¬ 
tion of FPMB to be 3% greater than the normalization 
of the FPMA. For NGC 1320 the FPMB normalization 
is 9% lower. Balokovic et al. (2014) also investigate the 
cross normalization between the two FPMs for NGC 424, 
NGC 1320 and IC 2560, and our results agree. Full de¬ 
tails of NuSTAR in-orbit calibration is presented in Mad¬ 
sen et al. (in preparation). 


4. SPECTRAL FITTING RESULTS 

Our spectral fitting with the TORUS model reproduces 
the data well in all ten sources, with reasonable reduced 
X 2 . We confirm the Compton-thick nature of all 10 AGN, 
with Ah constraints in excess of 1.5 x 10 24 cm -2 . The 
Compton-thick nature of 2MFGC 2280 is shown here for 
the first time. 

The best-fit spectral parameters are presented in Table 
[3j along with their 90% confidence intervals calculated 
using a Ax 2 =4.61 criterion on two interesting parame¬ 
ters in order to minimize degeneracies. The uncertainty 
on the luminosity is propagated from the uncertainty on 
the normalization of the TORUS model. Details of the 
parameters of the components used to fit the soft X-ray 
excesses are listed in Table |4] 

The best fit 0 t or measured in our sample spans a large 
range, from 26-80°, limited by the allowed range of the 
model. The uncertainties on d t0 i range from 5-40°, with 
the best constraints coming from NGC 1068, NGC 4945 
and Circinus which have the highest signal-to-noise spec¬ 
tra in the sample. While the fit to Mrk 34 gives a fairly 
large 9 t0 r, the constraints are poor, with statistically al¬ 
lowed values ranging from 26-80°. The case is similar for 
2MFGC 2280, albeit with slightly tighter constraints. 

In two sources, NGC 3079 and NGC 4945, we discover 
bimodality in statistic space for 9 toI1 with local minima 
at small and large opening angles. We show their x 2 as 
a function of d t0 r in Figure [5} For NGC 3079, the large 
opening angle provides a marginally better fit, whereas 
for NGC 4945 the small opening angle is a significantly 
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Figure 5. y 2 as a function of 4.or for fits to NGC 3079 and 
NGC 4945, showing local minima at both small and large open¬ 
ing angles. For NGC 3079 the large opening angle is marginally 
favoured, while in NGC 4945 the small opening angle is signifi¬ 
cantly favoured. The dashed lines indicate the minimum (26°) and 
maximum (84°) ©tor values allowed by the model and the dotted 
lines show the 68, 90 and 99% confidence levels which correspond 
to Ay 2 values of 1.00, 2.71 and 6.63 respectively. 


better fit. 

For two sources, NGC 1386 and NGC 4945, we find 
that the best-fit 0 to r reaches the lower limit on this pa¬ 
rameter (26°), which implies that both the sources are 
highly covered. We also find that the T values for these 
fits deviate significantly from the canonical value of 1.9 
(steep at 2.52 for NGC 1386 and flat for NGC 4945 at 
1.58). The fact that the deviations from the canonical 
value occur at the edge of the 9 tor parameter space sug¬ 
gests that the true 9 tor of these sources lies outside of the 
r ange of the TORUS m odel (i.e. d t or<26°). 


Brightman & Nandra (2011a) also present a model of 
a fully covered source (i.e. 0° opening angle), which in¬ 
cludes variable iron and elemental abundances, known 
as sphere. We fit this model, with the iron abundance 
free, in place of the TORUS model for NGC 1386 and 
NGC 4945, and present the results in Table [3] This is 
also done for NGC 3079 since the small opening angle 
also provides a good fit to that source. An improve¬ 
ment in the fit statistic is found for both NGC 3079 and 
NGC 4945, where the spherical model produces a sig¬ 
nificantly better (>90% significance) fit than the TORUS 
model for both (Ax 2 =4 and 136, respectively). The iron 
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Table 2 

Details of the cross-normalizations. 


Source name 

FPMB/FPMA 

XMM-Newton/ FPMA 

Models used 

(i) 

(2) 

(3) 

(4) 

NGC 424 

1 07+ u - i;i 

A * u ' -0.12 

0.971";- 

TORUS+POWERLAW 

NGC 1068 

1 03 +0 03 
±.uo_o 02 

- 

TORUS+POWERLAW+ZGAUSS+ZGAUSS 

2MFGC2280 

1 02 +0 ' 16 
1 - uz -0.14 

- 

TORUS 

NGC 1320 

°- 91 ±0.07 

0-89±°;l® 

TORUS+POWERLAW+ZGAUSS 

NGC 1386 

o-so±S:?S 

°-84±o.i6 

TORUS+POWERLAW+APEC+ZGAUSS 

NGC 3079 

1 06 +011 
A * UD -0.10 

0-65i°;i! 

TORUS+POWERLAW+APEC 

IC 2560 

1 07 +0 ’ 09 

A,u ' —0.09 

0-82±°;l«, 0.92±S;g 

TORUS+APEC+ZGAUSS 

Mrk 34 

1 09 + °' 27 
A * uy -0.22 

0.93±°;i?, 0.8il“ : 4 I 

TORUS+POWERLAW+APEC 

NGC 4945 

1 02 +0 ' 02 
l.UZ—o 02 

- 

TORUS+POWERLAW+ZGAUSS 

Circinus 

1 03 + °' 01 
i.uo-o oi 

- 

TORUS+POWERLAW+ZGAUSS 


Note. — Column (2) gives the cross-normalization between the two FPMs, column (3) gives the cross-normalization between FMPA 
and XMM-Newton (MOS 1 and MOS 2 respectively for IC 2560 and pn and combined MOS for Mrk 34) and column (4) lists the models 
used to fit each of the spectra. 


abundance has been estimated to be slightly higher, al¬ 
beit statistically consistent with solar metalicity. The T 
inferred from the sphere model is also consistent with 
the canonical value for both sources. For NGC 1386, the 
SPHERE model produces extreme T and iron abundances 
values and iherefore it is not likely to be a good descrip¬ 
tion of the spectrum. 

Due to the range of 8 tor allowed for the TORUS model, 
we cannot test if a very large opening angle (^90°) would 
provide as equally good a solution for NGC 3079 or 
NGC 4945. However for the local minimum at the largest 
angle, the T for NGC 4945 is even flatter (1.43) than the 
smallest opening angle and does not suggest that a large 
opening angle would describe the spectrum well. For 
this solution, Vh=2.64x 10 24 cm -2 and Lx=2.27xl0 42 
ergs -1 . 

We conclude from this that NGC 3079 and NGC 4945 
are heavily buried sources with a covering factors close to 
unity. For all further analysis we assign a value of 0 tor =O° 
to these sources. We discuss this result later, especially 
for NGC 4945 where previo us results have favou red a 
very low covering factor (e.g. Madejski et al.pOOO ). 

The best-fit unfolded spectra tit with the TORUS model 
are presented in Figure [6j This illustrates the variety of 
spectral shapes seen above 10 keV for these Compton- 
thick sources, as well as the range in signal-to-noise in the 
sample. The data-to-model ratios are presented in Figure 
0which indicate how well the TORUS model fits the data. 
The worst reduced \ 2 from the TORUS fits result from the 
spectrum of NGC 1068 (1-72), where there is significant 
curvature in the data-to-model ratio, implying that the 
TORUS model is not a good description of the data. This 
was also the conc lusion from the detailed analysis of this 
source in Bauer et al. (20141, which utilizes data from a 


large array of telescopes and several model combinations. 
They find that the data require a complex combination 
of models to reduce this curvature and that a monolithic 
torus structure is not likely. 

First in our investigation into what drives the cover¬ 
ing factor in AGN, we examine how the measurement of 
$tor in our analysis relates to the other measurements 
made with the TORUS model, namely Nr and T, by 
plotting these quantities a gai nst each other along with 
their uncertainties (Figure [8]). No uncertainties can be 


derived for the f?tor=0° values assigned to NGC 3079 
and NGC 4945 as these have been determined from the 
SPHERE model which has a fixed opening angle of zero. 
There is no clear relationship between these quantities. 

5. POTENTIAL BIASES AND SYSTEMATICS 

In our analysis, we fix the torus inclination angle to an 
edge-on position of 87°. This is primarily motivated by 
the observation that the inclination angle does not have 
a large effect on the observed spectrum above 10 keV and 
thus allows us to reduce the number of free parameters in 
the fit. Furthermore an edge-on inclination angle allows 
the full range of opening angles to be explored. Fixing 
the inclination angle also avoids the scenario where the 
inclination angle approaches the opening angle, which 
produces a partially covered solution due to the angular 
binning in the TORUS model. We investigate whether 
allowing the inclination angle to be free allows it to be 
constrained. This is n ot the case , ho wever Circinus is an 
exception, as found by Arevalo et al. (|2014|), thus we have 
allowed the inclination angle to be free for this source. 

The NuSTAR band is also well suited to studying the 
high-energy cut-off in the X-ray spectra of AGN, and has 
been used to constrain this in a number of bright unob¬ 
scured sources (e.g. SWIFT J2127.4+5654, IC 4329A, 
3C 382, M CG-5-23-16, |Marinucci et al.12014 Brenneman 


|et al.|[2014[ |Ballantyne et al.||2014 Bafokovid et al. sub¬ 

mitted, respectively). However, the TORUS model does 
not include this parameter and thus we have neglected its 
effect in this work. The potential effect of neglecting the 
high-energy cut-off here is to over-estimate the Comp¬ 
ton hump, which is produced by the down-scattering of 
high-energy photons. As the covering factor roughly cor¬ 
relates with the strength of the Compton hump, then 
over-estimating the Compton hump will lead to a sys¬ 
tematic under-estimation of the covering factor. We ex¬ 
pect this effect to be very small as the lowest high-energy 
cut-off energy detected so far by NuSTAR is at 1 08 keV 
in SWIFT J2127.4+5654 (|Marinucci et al.||2014|), which 
is far above that of the Compton hump 30 keV). For 
photons originating around a cut-off of this energy to af¬ 
fect the Compton hump, at least 10 scatterings would 
be required. To understand this effect fully, the TORUS 
model should be improved to include the high-energy cut- 
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Energy (keV) 


Energy (keV) 


Figure 6. Unfolded spectra for all sources fit with the TORUS model, where NuSTAR FPMA data are shown in red, FPMB data are shown 
in green and XMM-Newton data are shown in black. The solid black lines represent the sum of all model components, while the dashed 
line shows the TORUS model. Dot-dashed lines represent the secondary power-law model and the dotted lines represent the thermal plasma 
model APEC, both used to fit the soft-X-rays. The spectra have been binned for plotting purposes, so that each bin has a detection of at 
least 3-cr up to a maximum of 10 data points per bin. 
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Table 3 

Best-fit spectral parameters of the TORUS and sphere models. 







TORUS fits 





Source name 

redshift 

PHA bins 

x 2 

x 2 

A h 

r 

$tor 

logio^x 

logicRx 

(i) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

NGC 424 

NGC 1068 

2MFGC 2280 

0.0118 

0.0038 

0.0150 

121 

1140 

36 

139.27 

1950.47 

13.18 

1.23 

1.72 

0.43 

5.41++ 7 

6 32 +0 ' 44 
u -° -0.05 

2 ‘50+ 0,58 
z,ou - 1.00 

2 20 +U-34 
—0.16 
901 +0.04 
z,ol -0.05 

1 96” 1 " 0,26 

78-19+5+ 

58.99+75 

78.97++ 

-10.80 

-10.51 

-10.98 

43-96++ 
42.87+06 
43 26+ 0,31 

‘TO.ZO— q 74 

NGC 1320 

NGC 1386 

NGC 3079 

0.0089 

0.0029 

0.0037 

172 

71 

214 

190.37 

77.46 

270.38 

1.18 

1.31 

1.33 

ioo.oo+ 0 . 81 

5 61 + 2,11 

9 07 + 0.36 
z -0.25 

1 cc+ 0.21 
L66 -0.19 

9 qq + 0.08 

z,yz -0.44 

1 69 +0,11 

1 - Dy — 0.15 

60 17+ 10 - 69 
DU,i ‘ _ 14 . 91 
99 cq+10.79 
oo.<jo_7.23 

79 75+ 1,73 

/y. /o_9 36 

-10.95 

-11.68 

-10.78 

42 79+ 0,12 
* y -0.09 

41 84+ 0,26 

4i,O4 -0.05 

41 96+ 0,19 
4i.» 0 _o. 33 

IC 2560 

Mrk 34 

NGC 4945 

Circinus 

0.0096 

0.0510 

0.0019 

0.0014 

165 

69 

906 

1715 

182.51 

71.03 

1108.75 

1880.27 

1.19 

1.22 

1.24 

1.10 

100.00tSg.84 

50 43+ 49 - 57 
ou.^o_ 31 01 

2 54+ 0,13 
z -°^-0.l5 

4.80-0.42 

9 kq+ 0.20 
z,oo -0.20 

1 70 +I .25 
—0.56 

1 58+ 0 - 04 

i - OO -0.04 

9 97 +O.O 5 
z,z ' —0.07 

59 08+ 9 ' 57 

oU.UO — 20 Q5 

79 n +6.89 
' z - 0i — 46. 51 
26.00+ 50 

33.79±i;|| 

- 11.66 

-11.52 

-9.65 

-9.63 

42 95 + 0,11 

4Z.UO-0 13 

44 18+ 0,39 

4^.io_ 0 .3 8 

41 92+ 007 

0 os 

42 51 +0.07 

TZ.Ol — 0.09 






SPHERE fits 





Source name 

redshift 

PHA bins 

X 2 

x 2 

A h 

r 

Fe abund. 

logio^x 

logio^x 

(i) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

NGC 1386 

0.0029 

71 

73.74 

1.25 

0 76~*~ 9 ' 45 

U. /O-0.19 

3.00+ 

10 . 00+ 44 

-11.56 

39.67++ 

NGC 3079 

NGC 4945 

0.0037 

0.0019 

214 

905 

266.42 

972.64 

1.31 

1.09 

1 84”^°' 32 

9 o ^ - ! -9 - 24 
z,zo -0.07 

J-. 8 D_o 23 

1 7Q+0.13 

1 •— 0.06 

1 42+ 0 - 34 
- L,ftz - 0.43 

1 qo+0-16 
l.oo-o 37 

-10.83 

-9.66 

41 53 + 0,45 

41.00 — 0 43 

42 05+°' 17 
4Z.UO_o 1Q 


Note. — Parameters determined by the TORUS model are listed in the top rows. For NGC 1386, IC 3079 and NGC 4945 we also present 
the best-fit parameters from the sphere model as the TORUS model finds that these sources are highly covered (bottom rows). Column (1) 
lists the source name, column (2) gives the redshift of the source, column (3) shows the total number of pulse height analysis (PHA) bins 
used in the spectrum, column (4) gives the x 2 of the fit, column (5) gives the reduced x 2 , equal to x 2 divided by the number of degrees of 
freedom in the fit, column ( 6 ) gives the Nn in units of 10 24 cm -2 , with uncertainties, column (7) gives the photon-index, column ( 8 ) gives 
0t or , determined by the TORUS model in degrees or the iron abundance with respect to solar hydrogen abundance from the sphere model, 
column (9) gives the logarithm of the observed 3-79 keV NuSTAR flux in erg cm -2 s —1 and column (10) gives the logarithm of the intrinsic 
(deabsorbed) rest-frame 2-10 keV luminosity in ergs -1 . +u indicates that a parameter has reached the upper limit when estimating the 
uncertainty. —I indicates that the lower limit has been reached. 


Table 4 

Best-fit spectral parameters of the secondary power-law, APEC and Gaussian line models. 


Source name 

r 

Api/10 - 5 

Tapec 

Aapec/10 - 5 

^line 

Aline/10 -5 

(i) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 




TORUS fits 




NGC 424 

3 7 Q + 3 45 

0 . < y_ 0 .93 

11.27+;?! 

- 

- 

- 

- 

NGC 1068 

2 29 +0 ‘ 10 

z - zy -0.09 

163.20+?+ 


- 


- 

2MFGC 2280 


- 


- 


- 

NGC 1320 

3 70+°- 46 
°- /u -0.41 

12 91+ 3 - 11 
iz.yi _2 51 


- 

6 55 ”*" 9,95 
°'0°—o.04 

0 57 + 9 - 29 

U * 0 ' -0.17 

NGC 1386 

2 56+ 1 ’ 72 
z.oo _2 76 

1 61+ 1 ' 07 
i-oi-l.06 

n hr 4-9-98 

U ‘ b °-0.09 

7 01 +1.06 

' - 01 —3.57 

6 42+ 9 - 04 
u, ^ z -0.06 

0 47 +°- 17 

u ‘^ * —0.17 

NGC 3079 

1 95 + 1 ’ 70 

i.yo _3 34 

0 qo+9.17 
o.oo _3 38 

0.93+;?! 

9 92 - * -21,45 
-5.72 


- 

IC 2560 

- 

- 

0 . 68 +;?? 

9 90 +I.I 6 
Z,ZO -0.62 

6 43 + 9 - 02 
o -^°- 0 . 0 i 

0-67±°;ig 

Mrk 34 

2 90 + °' 81 
z - yu -0.28 

1 q 8 +°- 99 

i.uo-l 19 

0.82+;?! 

0 92 +0,49 

U - yz -0.54 

- 

- 

NGC 4945 

i-90+;f 8 

90 79+38.76 

z y .' Z — 16.28 

- 

- 

- 

- 

Circinus 

2 . 05 +:?! 

226.28+?? 8 2 ? 7 - 

- 

6 57 + 0 04 
D *°' - 0.00 

12.86li;i5 




SPHERE fits 




NGC 1386 

1 20+ 0 ’ 33 
1 * zu —0.30 

1 70 +O .85 
i *'°-0.63 

0 . 68 +;?! 

7 60 +1,58 
'-00-1.52 

6 39 + 0 - 05 

u.oy _ 0 04 

0 . 54 +;?? 

NGC 3079 

1 qo+0.37 

J--o 8 _o 25 

6.09+;!? 

0.89+;?? 

11-82+;?? 

- 

- 

NGC 4945 

1 OO+ 022 
l.UU-o 25 

12.57+;?? 

- 

- 

- 

- 


Note. — Column (1) gives the source name, column (2) gives the power-law index of the secondary power-law, column (3) gives the 
normalization of the secondary power-law, column (4) gives the temperature of the APEC model in keV, column (5) gives the normalization 
of the APEC model, column ( 6 ) gives the energy of the Gaussian line in keV and column (7) gives the normalisation of the Gaussian line. 
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Figure 7. Data-to-model ratios of the fits with the TORUS model. 
The colors are the same as Fig. 5. The binning for this figure has 
been increased in order to better see deviations and trends in the 
residuals. 


uncertainties for all of these except NGC 424. For this 
source, XMM-Newton measures Ah= d.hlJ'yXlO 24 cm -2 
and NuSTAR measures Ah= 2.Up® x 10 24 cm -2 . Allow¬ 
ing for the Ah to vary like this does not however change 
our result on 9 tor . 

Lastly we investigate how further spectral components 
added to the fit may affect our results. Specifically, we 
investigate the addition of the semi-infinite slab reflec¬ 
tion model pexrav in order to represent a contribution 
to reflection from the accretion disc. We use the pure 
reflection component of this model and fix T to that of 
the TORUS model, set the high-energy cut-off to max¬ 
imum (1 x 10 6 eV), the adundances to solar and the 
cosine of the inclination angle to 0.45. The normaliza¬ 
tion is free. In four sources, NGC 1068, 2MFGC2280, 
NGC 1320, NGC 1386 and IC 2560, the normalization 
of the PEXRAV component falls to zero in the fit. For 
NGC 3079 and NGC 4945, the addition of the PEXRAV 
component has no effect on (9 tor , whereas for Circinus, 
the change is —5°. For NGC 424 the addition of pexrav 
to the fit causes a shift in 9 t or of —25°. While the change 
in 0 t or is large for NGC 424, when considering the sample 
as a whole we conclude that the addition of this compo¬ 
nents does not appear to introduce a systematic effect 
in the determination of 9 tor . We note that the PEXRAV 
component added here is unabsorbed and not subjected 
to the same absorption as the primary power-law. To 
treat this correctly, the PEXRAV spectrum would need to 
be added to the intrinsic emission of the Monte-Carlo 
models that calculate the effect of Compton scattering 
within the TORUS model, which is not easily done. 


off. 

For this study, we have a heterogeneous sample in 
terms of X-ray spectral coverage and most of our sample 
have accompanying soft X-ray data to constrain the spec¬ 
tral components below 10 keV. However, for four sources, 
we use NuSTAR data alone. No soft X-ray data exist 
for 2MFGC 2280, meanwhile for NGC 1068, NGC 4945 
and Circinus, the NuSTAR data are high signal-to-noise 
enough that they can constrain the TORUS parameters 
alone. We test this assumption for all three sources by 
adding in an XMM-Newton observation to model the soft 
X-rays below 3 keV. We find an additional APEC model 
is required to fit these data in addition to the TORUS 
and power-law model already in place. Fitting the data 
simultaneously with the soft component constrained by 
XMM-Newton produces only small changes in 9 tor , and 
within measurement uncertainties. We conclude that the 
exclusion of the soft X-ray data does not bias our results 
significantly in these sources. 

When fitting XMM-Newton and NuSTAR data to¬ 
gether, we have kept the spectral parameters fixed be¬ 
tween the two data sets. Since the observations were 
not made simultaneously, spectral variability could be 
missed, with the exception of the intrinsic luminosities, 
which is accounted for by the cross-normalisation be¬ 
tween the two observatories being a free parameter. We 
investigated spectral variability in the five sources where 
we use both XMM-Newton and NuSTAR data together 
by allowing the Ah parameter to vary between data sets. 
We find that the Ah parameter is consistent within the 


6 . COMPARISON TO PREVIOUS RESULTS 


The use of the TORUS model here to determine the 
opening angle of the torus is fairly novel, as is the use 
of X-ray torus models to understand high-energy emis¬ 
sion in heavily obscured AGN. Therefore, in Section [2] we 
compared the spectral fit parameters determined from 
the classical PEXRAV model and the more recent MY- 
TORUS model when fitted to simulated spectra from the 
TORUS model for a range of 9 tor and Ah- Here, we 
compare the fits to real NuSTAR data with the TORUS 
model presented here to those made with PEXRAV and 
MYTORUS from previous published works. We list the 
different Ah and T values obtained from each in Table [5| 
In our initial analysis, we found that PEXRAV system¬ 
atically underestimates T with respect to both TORUS 
and MYTORUS, most likely because it cannot reproduce 
the full strength of the Compton hump produced by a 
torus geometry. This is also found in our comparison of 
fits to real data, as fits with PEXRAV produce a lower T 
than both the TORUS or MYTORUS model for NGC 424, 


NGC 1068, NGC 1320 and IC 2560. 

As for fits with MYTORUS with respect to the TORUS 
model, the T measurements are consistent with each 
other within the uncertainties, with the exceptions of 


NGC 1068 and NGC 4945. For NGC 1068, Bauer et al. 


(20141 use a large, combined set of X-ray spectral data 
from different observatories with far more detailed spec¬ 
tral modelling than is done here, including modelling of 
the host spectrum. Futhermore, the MYTORUS fits are 
done in decoupled mode. For NGC 4945, the TORUS 
model measures T ~0.2 lower than MYTORUS. However, 
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Figure 8. Best-fit measurements of #tor against best-fit measurements of TVh and F from the TORUS model. NGC 3079 and NGC 4945 
are best-fit by the sphere model implying that these sources have 0 tor =O°. No uncertainties can be calculated for this, however, as the 
SPHERE model has a fixed opening angle of zero. 



the m easurement of T using MYTORUS in Puccetti et al. 
(2014) was done on a time-resolved basis with MYTORUS 
in decoupled mode, where the parameters of the reflec¬ 
tion components are not fixed to those of the transmit¬ 
ted component. However, we do measure a lower TVh 
for NGC 4945 than the MYTORUS model, which may 
be the reason for the discrepancy. Furthermore, a fit 
with the SPHERE model produces a better fit to the data 
and a T value consistent with MYTORUS (1.78+o)o6)- For 
Mrk 34, the difference in the measured T values between 
the two models is ~ 0.5, with MYTORUS producing the 
steeper slope. The TORUS model, however, measures 
Ah> 10 25 cm -2 in Mrk 34, whereas MYTORUS is lim¬ 
ited to Ah< 10 25 cm -2 , and, indeed the upper limit is 
reached in the fit. The inability of MYTORUS to mea¬ 
sure columns greater than this is the likel y cause o f dis¬ 
agreem ent in Mrk 34, as also noted by Gandhi et al. 
(2014). The two measured values are nevertheless con- 


smeared out by reflection and that a low covering factor 
would be requ i red to produce such weak reflection. 


sistent within the large uncertainties in T. 

6.1. NGC 4945 

The high covering factor determined here for 
NGC 4945 is in disagreement with previous analyses 
which conclude t hat NGC 4945 has a small covering fac¬ 
tor, of order 0.1 ([ Madejski et al.||2000[ 


Yaqoob 2012 F 
drawn from th 


1 


Done et al. 


2003 


rom the tact that 4945 is variable above 

10 keV and that it has a relat ively weak reflected com¬ 
ponent. Madejski et al. (2000) used Monte-Carlo simu¬ 
lations to show that the fraction of unscattered photons 
reaching the observer is much higher for low covering fac¬ 
tors (63% for 0 tor =8O°) than high covering factors (19% 
for 0 t or=lO°). The argument states that with an optical 
depth to Compton scattering of 2-3, the intrinsic vari¬ 
ability of the AGN with a high covering factor would be 


Done et al. (2003) also conclude that the Compton- 
thick material cannot cover the whole source in 
NGC 4945 as the 6.4 keV emission is spatially extended 
in Chandra observations and m ust be illuminated by hard 
X-rays from the AGN. Done et al. (2003) also suggest 
that NGC 4945 must be fully covered due to the lack of 
high ionization optical/IR lines, although this material 
does not n eed to be Compton-thick. 


Yaqoob (2012) presented an extensive broad-band X- 
ray spectral analysis of NGC 4945 using Suzaku , Bep- 
poSAX and Swift /BAT data using MYTORUS, TORUS and 
SPHERE models and also prefer the small covering factor 
solution due to the hard X-ray variability. 

A detailed spectral and tempor al analysis of NGC 4945 
with NuSTAR was presented in Puccetti et al. (20141 
utilizing the MYTORUS model in decoupled mode. The 
hard X-ray variability was confirmed, where variations 
of a factor of two above 10 keV were reported. They 
estimate the covering factor for this source, using the 
ratio of the reflected component to the direct component. 
Since this is very small, they conclude that the covering 
factor is ^0.13, and fairly constant with flux. 

Using the TORUS model, we have found that the X-ray 
spectral shape of NGC 4945 can also be produced by a 
source with a high covering factor of obscuring Comp¬ 
ton thick material and that for this model this in fact 
provides a better fit to the X-ray spectrum than the low 
covering factor scenario. Figured] shows y 2 as a function 
of 0tor for the fit to NGC 4945. Although this shows 
a local minimum at large opening angles (low covering 
factor), the minimum at small opening angles (large cov¬ 
ering factor) is significant ly lower. Indee d , with their 
Monte-Carlo simulations, Madejski et al. ( 2000|) using 
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Table 5 

Comparison of fits with the TORUS model to previous analyses with pexrav and MYTORUS 


Source name 

N h (mytorus) 

Nh (torus) 

T (pexrav) 

T (mytorus) 

r (torus) 

$tor 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

NGC 424 

NGC 1068 

NGC 1320 

IC 2560 

Mrk 34 

NGC 4945 

3±1 

10±4.4 

4i^ 

10 ±“ 

9 4A+' U 
z - 40 —1.08 

3.5±0.1 

5 - 41 + 123 y 

6 32 +0 ' 44 

°- oz -0.05 

ioo.ool“ 081 

100.00lSg.84 

50.43l«; s 0 ( 

2-54l°;i| 

1.66T0.09 

1.57T0.02 

1.3T0.1 

9 9+O.I 
Z,Z -0.2 

2 07+O- 11 

z - u/ -0.09 

9 9p»+0.02 
z,zu -0.01 

1 . 6 d= 0.2 

2.55 (f) 

9 9+O.2 
z,z —0.3 

1.77-1.96 

2 20”*" u ' a4 

z.zu—0.16 

9 qi +0.04 
z ’ o± -0.05 

-i aa +0.21 

9 cq+0.20 

Z.OO — 0.20 

1 70 + 1.25 
1-/O-0.56 
i co+0.04 
± - oo -0.04 

78.19l)’ 5 4 g 1 

58.9912Y6 

eo.utiK'®) 1 

59.08l| 0 5 J 5 

79 n +6.89 
iz.oi_ 46 51 

26.00+ 50 

Circinus 

10 .Oil .8 

4.85l°; 3 4 9 

- 

2.34±0.02 

9 97+O.O5 
z - z ‘-0.07 

33.79li'H 


Note . — Frompa lokovic et al. (2014) for NG C 424, NGC 1320 an d IC 2560,|Bauer et al.| ^2014h for NGC 1068,|Gandhi et al.|(|2014|) for 
Mrk 34, |Puccetti eT al. (2U14) tor Ngg 4945 and | Arevalo et ah] ([2014 ) for Circinus. Uolumn (l) lists the source name, column [ 2 ) gives the 
Nn measured by MYTORUS in units of 10 24 cm -2 . +u indicates that the upper constraint on the Nn from MYTORUS is beyond the upper 
limit of the model. Column (3) gives the Nn measured by the TORUS model in the same units, column (4) gives the photon index measured 
by pexrav, column (5) gives the photon-index measured by MYTORUS, where (f) indicates this parameter has been fixed, and column ( 6 ) 
gives the photon index measured by TORUS. Column (7) lists 0tor measured by the TORUS model in units of degrees. -\-u indicates that a 
parameter has reached the upper limit when estimating the uncertainty. —I indicates that the lower limit has been reached. 


RXTE data also find that the spectral shape favours a 
small opening angle, finding y 2 =68.5 for #tor=20° and 
y 2 =75.4 for 0 tor =8O°. Their conclusion regarding the 
large opening angle is instead driven by the variability 
and their Monte-Carlo simulations. However, it is not 
clear whether these results based on the variability are 
dependent on the geometry of the torus used (a torus 
with a square cross-section), or if these results are en¬ 
ergy dependent. 

Our results also imply an optical depth to scattering of 
1.5-1.7, which is lower than previously considered, and 
thus the effect of scattering is slightly diminished. How¬ 
ever the fraction of directly transmitted photons is still 
low when considering a high covering factor, which is 
hard to reconcile with the hard X-ray variability. 

The conclusions regarding the covering factor of 
Compton-thick material surrounding NGC 4945 drawn 
from the variability and those drawn from the spectral 
shape are at odds and the models used to draw these 
conclusions both have their limitations. Concerning the 
variability, it has been assumed that scattered photons 
cannot transmit the intrinsic variability of the source. 
However, Compton scattering prefers forward (and back¬ 
ward) scattering with small angles, so the difference in 
light travel time between transmitted and scattered pho¬ 
tons need not be large. The difference in light travel 
time is also dependent on the distance of the scatterer 
from the central source. For the TORUS model, the anal¬ 
ysis is limited to a range of covering factors, not allowing 
investigations of very small covering factors (~ 0.1) or 
very high covering factors (~ 0.9 — 0.99), and also does 
not allow one to decouple the transmitted and scattered 
components. The assumption of a smooth matter distri¬ 
bution is also unlikel y to be accurate as it is most likely 
to be clumpy (|Yaqoob|2012). It is clear that further work 
is required to fully"understand the nature of the absorber 
in NGC 4945. 

7. WHAT DETERMINES THE COVERING FACTOR OF 
THE OBSCURER IN AGN? 

The obscurer in local AGN is widely regarded to be 
a cold molecular torus, which man y recent resu l ts im¬ 
ply has a clumpy constituen cy (e.g. Elitzur & Shlosman 
2006 Markowitz et al.||2014). As discussed, high-energy 


X-rays are ideally suited to the study of the obscuring 
material, since Compton scattering effects from the gas 
in the obscuring medium dominate in this regime (>10 
keV). Our fits to NuSTAR data with the TORUS model, 
which assumes a smooth torus, support the general torus 
paradigm. In our analysis we determine a wide range of 
9 tor allowed by the model, and furthermore identify three 
sources where the spectral fits indicate small or zero 9t or - 

Although our sample is small, we investigate what 
could be physically influencing the opening angle. As 
the obscured fraction is known to depend on X-ray lu¬ 
minosity, we first investigate how the covering factors 
derived here depend on Lx, using the 2-10 keV band 
and the intrinsic luminosity determined from the model. 
Figure [9] plots these quantities against each other with 
their measured uncertainties. A strong anti-correlation 
is seen, as expected, where one of the most luminous 
sources, NGC 424 with logi 0 (Tx/ergs _1 )=43.96lU+ 
has a small, relatively well constrained covering factor 
(/ c =0.20±g;g®), while Circinus with a moderate luminos¬ 
ity of logio(L x / erg s _1 )=42.5llH;Hg has a larger covering 
factor of / c =0.83to;o2- Such a correlation is expected due 
to more luminous sources sublimating dust in the inner 
edge of the torus at larger distances. For the same verti¬ 
cal extension of the torus, a larger radius of the inner part 
of the torus gives a lower covering factor of the central 
source due to geometrical effects. 

We fit a simple linear model, y = mx+c to the covering 
factor vs. logio^x data using the IDL function LINFIT, 
which utilises y 2 minimisation and takes into account the 
uncertainties in the covering factor. We find that f c = 
(—0.41±0.13)logio(Tx/ erg s _1 )+18.31± 5.33, where the 
uncertainties are 1-a. We plot this function along with 
the uncertainties in Figure]!)} We do not include the data 
for NGC 3079 or NGC 4945 in the fit as these have no 
uncertainties in the covering factors assigned to them. 
Notably, however, they both agree very well with the 
fitted model. 

We also compare the covering factors to the obscured 
fraction of local AGN from three studies. We compare 
to rec ent o bscured fractions present ed by |Burlon et al.| 
poTTf and |Vasudevan et al.j (|2013|), both har d X-ray 
selected samples from Swift/BA1', and that of Bright | 
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Figure 9. The covering factor as a function of intrinsic 2-10 keV luminosity. Lv. derived in our analysis, where the best-fit linear model 

f c = (—0.41 =t 0.13)logio(Lx/ergs~ x ) + 18.31 ± 5.33 is plotted. NGC 3079 and NGC 4945 arc excluded from the linear function fit since 

they both have no uncertainties on their covering factor. Nonetheless, they both agree well with the derived function. 


man & Nandra| (|2011b|, a mid-infrared, IRAS -selected 
sample. Ail three obscured fractions are defined as the 
fraction of sources with Ap> 10 22 cm -2 , calculated in 
different luminosity bins. We plot these obscured frac¬ 


tions in Figure 1 10[ The uncertainties in the Brightman 
& Nandra ( 2011 bp data points are binomial. The Burlon 


et al. (|2U11 ) curve is calculated by dividing the X-ray lu¬ 


minosity function (XLF) of obscured A GN by the total 
XLF, done in the 15-55 keV band. The Vasudcv an et al.| 
(2013) line is a running average using 30 sources per bin. 
The two different lines for this sample are derived from 
where there are uncertainties on the Ah measurement, 
and the upper Ah bound is used for the upper line and 
the lower Ah bound is used for the lower line. 

The obscured fraction determined in these three stud¬ 
ies agree very well with each other despite the differing 
selections and determinations, declining from a peak at 
Lx= 10 42-43 ergs -1 towards higher luminosities. All 
three studies also find evidence for a decline in the ob¬ 
scured fraction towards lower luminosities. 

Also in Figure [lOj we overplot the covering factors de¬ 
rived from our sample. We find good agreement between 
our derived covering factors and the obscured fraction 
for Lx> 10 42 5 ergs -1 . However, while previous studies 
have found evidence for a decrease in the covering fac¬ 
tor at low luminosities, we find that our sources in the 
Lx= 10 41-42 ergs -1 range are heavily buried in mate¬ 
rial with high covering factors. A larger more complete 
sample is required to show if this disagreement is statis¬ 
tically significant or due to the low number statistics of 
our sample. 


We note that the covering factors derived here are 
those of the Compton-thick gas surrounding the AGN, 
since the TORUS model assumes a constant density torus 
with a constant Ah as a function of inclination angle. 
Therefore any additional covering by Compton-thin gas 
would not be recognised, in which case the covering fac¬ 
tors here could underestimate the total covering factor. 
Due to the agreement between the Compton-thick cov¬ 
ering factor and the obscured fraction (for Lx> 10 42 5 
ergs -1 ), this does not appear to be the case in our 
sample. This implies uniform torus covering factors for 
Compton-thin and Compton-thick AGN given the same 
Lx, above 10 42 5 ergs -1 . Best estimates of the local 
Compton-thick fraction put it at ~ 20% (Burlon et al. 


2011 Brightman & Nandra 2011bI, which is lower than 


the covering factors determined nere, which could sug¬ 
gest that a larger population of CT AGN exists in the 
local Universe. 

Lastly, we briefly investigate what other AGN param¬ 
eters may be involved in determining the covering factor 
of the obscurer. We have already shown that the X- 
ray luminosity, which traces the bolometric power of the 
AGN and is thus dependent on the mass accretion rate, 
plays an important role. We next explore if the mass 
of the black hole, Mbh, or the fraction of the Edding¬ 
ton luminosity, AEdd, physically influences the covering 
factor. These quantities are notoriously difficult to de¬ 
termine in obscured AGN, as virial mass estimates from 
optical broad lines are not accessible. In these cases, the 
velocity dispersion of the stars in the bulge is often used 
to estimate Mbh from the Mbh— o'* relation, although 
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Figure 10. The covering factor of the torus as a function of intrinsic 2-10 keV luminosity, Lx> derived in our analysis in comparison to 
the obscured fraction of local AGN determined from Burlon et al. (2011), Brightman & Nandra (2011) and Vasudevan et al. (2013). The 
covering factors agree well with the obscured fraction for Lx> 10 42 ' 5 ergs -1 . However, at low luminosities, we find no evidence for the 
decline in the covering factor seen in those studies. 


there is evidence for la rge scatter in this relationship, 
especially at low mass (Greene et al. 2010). For some 
sources, water megamasers can provide robust black hole 
mass measurements. We assemble this data from the lit¬ 
erature and list it in Table [6] Some AEdd estimates are 
also available. The covering factor of the torus relative 
to these quantities are plotted in Figure [Tl] 

This preliminary investigation seems to show that the 
highest black hole mass systems in our sample have the 
smallest covering factor, while the smallest black holes 
have the highest covering factors, however the relation¬ 
ship is not statistically significant and a far larger sample 
is required to confirm this trend and to break the degen¬ 
eracy with Lx- As for the covering factor as a function 
of AEdd) our data suggest that high covering factors are 
exhibited in both low (10 -3 ) and h igh (~ 0.3) Asdd sys¬ 
tems, in line with the conclusions of Draper & Ballantyne| 
(2010) who find that CTAGN are made up of a composite 
population of both high (>0.9) and low (<0.01) Edding¬ 
ton ratio systems. 


8 . CONCLUSIONS 


In t his paper we have used the X -ray torus mod¬ 
els of Brightman & Nandra (2011a) and data from 
NuSTAH and XMM-Newton to determine the cover¬ 
ing factor of the Compton-thick gas in ten local CT 
AGN, NGC 424, NGC 1068, 2MFGC 2280, NGC 1320, 
NGC 1386, IC 2560, Mrk 34, NGC 3079, NGC 4945 and 


Circinus. We have also assessed the differences between 
the TORUS model, PEXRAV and MYTORUS. We find: 


• The slab reflection model, PEXRAV, does not eas¬ 
ily reproduce the Compton hump shape produced 
by a torus geometry, under-producing it for large 
covering factors, and over-producing it for small 
ones, resulting in a systematic offset in the param¬ 
eters obtained. We therefore discourage use of this 
model in the fitting of high-energy X-ray emission 
of CTAGN. Our results compare well with MY¬ 
TORUS for 1Vh< 10 25 cm -2 , where that model is 
valid, however we support the use of the covering 
factor as a free parameter in torus models. 

• Measurements of 9 tor are in the range of 26-80°, 
limited by the allowed range of the model, with 
uncertainties on these measurements ranging from 
5-40°. These correspond to covering factors in the 
range of 0.2-0.9, with uncertainties ranging from 
0.05-0.6. 

• The covering factor is a strongly decreasing func¬ 
tion of intrinsic 2-10 keV luminosity; when fitted 
with a linear function, we find f c = (—0.41 ± 
0.13)logio(L x / ergs -1 )+18.31 ± 5.33. 

• The individual covering factors derived here agree 
well with the average covering factor of local AGN 
as measured by the obscured fraction as a function 
of Lx above 10 42 5 ergs -1 . However, while previ¬ 
ous studies have found evidence for a decrease in 
the covering factor at low luminosities, we find that 
our sources in the Lx= 10 41-42 ergs -1 range are 
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Table 6 

Accretion parameters of our sample 


Source name 

(i) 

class 

(2) 

fc 

(3) 

logioix 

(4) 

logloM B H 

(5) 

logloAEdd 

(6) 

ref 

(7) 

NGC 424 

Sylh 

0 - 20 - 0 .o! 

43.96lS:w 

7.78 

-1.30 

1 , 6 

NGC 1068 

Sylh 

n 52 +0 - 04 

u - oz -0.07 

42.871°;“ 

7.59 

-1.42 

2 , 2 

2MFGC 2280 

Sy2 

°-19±o.09 

43.26l°“ 

0.00 

-99.00 

0 , 0 

NGC 1320 

Sy2 

0.501°;^ 

42-791 

7.29 

-1.54 

2 , 2 

NGC 1386 

Syli 

0.831°;°® 

41.841°;^ 

7.42 

-2.92 

2 , 2 

NGC 3079 

Sy2 

i nn+°- 16 
i.UU — 0,03 

4i.96l°;i° 

6.30 

-2.68 

7, 2 

IC 2560 

Sy2 

0 51+ 026 

U.O1-0.15 

42-951°;“ 

6.45 

-0.50 

3, 6 

Mrk 34 

Sy2 

0 Qn + 0 - 60 
U.OU- 0.12 

44.l8l° 0 ; 3 3 ° 

7.90 

-99.00 

5, 0 

NGC 4945 

FSRS 

i on - * -0 00 

41 QO + O- 07 
41.UZ-0 08 

6.15 

-99.00 

4, 0 

Circinus 

Sylh 

n QQ+ 0-01 
U - 8 ci - 0.02 

42 51 +0.07 

4Z.Ol o 09 

6.18 

-0.70 

8 , 9 


Note. — Column (1) lists the source name, column (2) gives the AGN activity class based on the optical spectrum, where Sylh indicates 
a Seyfert 2 with a hidden broad line region revealed in polarised light, Syli indicates a Seyfert 2 with a broad line region visible in the 
infrared, Sy2 indicates a Seyfert 2 with no evidence for a hidden broad line region and FSRS means a flat spectrum radio source. Column 
(3) gives the covering factor of the obscuring material derived in this work, column (4) gives logio(Lx/ergs -1 ) also derived here, column 
(5) lists the black hole masses from the literature in logarithm of sola r ma s ses, c olumn (6) gives the Eddington ratio published in the 
literature and column (7) list the references for these data: 1. |Bian &; Gu (2007 ) (black hole mass from stellar velocity dispersion), 2. 
|Marinucci et jd. (2012) ( black hole mass from stellar velocity dispersion} 3. |lshihara et al. (2001) (black hole mass from water megamaser 
emission), 4. [G rcenhill et al. (1997} ( black hole mass from wa ter megamaser emission), 5. puet al. (2008} (black hole mass from [OIII] line 
wi dth), 6. Baiokovic ct al. (20l4||; V. |Kondratko et al.|(|2005|); 8. Greenhill et al. (2003} (black hole mass from water megamaser emission); 
9. | Arevalo et al.| ( |2U14| ). 
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Figure 11. Our derived covering factor of the obscuring material compared to the black hole masses (in solar masses) and the Eddington 
ratios from the literature. 
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heavily buried in material with high covering fac¬ 
tors. A larger more complete sample is required to 
show if this disagreement is statistically significant 
or due to the low number statistics of our sample. 

• We find a conflicting result on NGC 4945, where 
our spectral analysis implies a large, almost 100%, 
covering factor, whereas previous results have con¬ 
cluded that this source has a very low covering fac¬ 
tor due to flux variability above 10 keV. We con¬ 
clude that model limitations in both cases are the 
likely cause of the disagreement. 
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